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Abstract 


Tha  diffusion  .erdel  of  Hutohlnson  at  al  (1971)  for  smII  partlclaa  has  baan 
extended  to  larga  partlelaa  in  a  vartloal  plpa.  Allowance  has  baan  made 
for  tha  slippage  In  tha  vartloal  direction  and  for  tha  r-ct  that  successive 
partlole  dlsplaoaaants  ara  no  loncar  antiraly  lndapandsnt.  It  U  shown 
that  tha  aquations  oan  bs  simplified  to  give  a  aiapla  expression  for  the 
diffusion  coefficient,  and  that  this  expression  gives  values  wnlch  are  in 
good  agreement  with  the  values  obtained  using  a  full  2D  simulation  or  the 
partlole  trsjeatorlea.  A  or  iter  Ion  Is  given  for  determining  when  a 
partlole  may  be  considered  to  be  "large". 
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t.  iinjiopocTiog 

The  motion  of  solid  partlolss  or  liquid  droplsta  In  a  turbulsnt  (as  flow  is 
a  problem  union  has  aany  praotloal  applications,  «(  In  pneummtlo  convoying 
and  spray  drying.  It  la  oftan  naeaaaary  to  account  for  part  tela  dlaparalon 
dua  to  turbulence,  as  a  simple  Nautonlan  calculation  uili  usually  under* 
astlaata  the  degree  of  dispersion.  Hutchinson  at  al  (1971 )  reviewed  the 
available  information  and  presented  a  nodal  which  traits  the  particle 
motion  as  a  diffusive  process.  This  model  was  shown  to  give  good  predic¬ 
tions  for  small  partlolss  (up  to  110  microns)  but  it  was  suggested  that  the 
model  may  not  be  applicable  to  large  particles. 

The  present  worn  Investigates  the  behaviour  of  large  particles  (greater 
than  1  mm)  In  a  vertloal  pipe  and  shows  how  the  Hutchinson  model  may  be 
adapted  to  allow  for  slippage  between  the  particle*  and  the  mean  gas  flow. 
This  report  also  shows  that  for  large  partlolss  there  Is  a  considerable 
degree  of  correlation  between  successive  particle  displacements.  A  method 
Is  presented  for  calculating  the  magnitude  of  the  correlation  terms,  and 
this  results  In  a  simple  expression  for  the  diffusion  coefficient.  The 
predicted  values  are  cornered  with  those  obtained  from  •  2D  simulation  of 
the  particle  trajectories. 

2.  DlffUSIOH  MODgL 

2.1  Diffusion  coefficient  for  small  particles 

The  model  of  Hutchinson  et  al  considers  the  particle  motion  In  two  parts 

(I)  constant  velocity  parallel  to  the  axis  of  the  pipe  and 

(II)  a  radial  motion  consisting  of  a  large  number  or  displacements  due  to 
partlole-eddy  Interactions.  This  radial  motion  Is  treated  as  a  dirfusive 
process  and,  assuming  cylindrical  symmetry,  the  density  distribution  <v  f  r .  r  . 
Is  glvsn  by  ths  solution  of  ths  diffusion  squstlon 

KV'W(r.t)  -  W(r,t)  •  3<r,t) 

where  S(r,t)  Is  the  source  term  and  n  Is  the  diffusion  coefficient  given  ry 

<v> 

K  -  v 

a 

where  <1  Is  the  meen  square  displacement  per  interaction  and  .  in  tie 
P 

Interaction  frequency.  This  expression  for  K  assumes  tnat  tne  ii'-*--!!  ns 
of  lndlvlduel  displacements  are  uncorrelateo  le  at  any  given  ••••  ;  c 

tide  has  an  equal  probability  of  moving  In  any  direction. 


t  - 


Asauaing  the  turbulence  la  hoMgeneous  ta  hat  no  preferred  direction,  man 

<«.  *>  •  2  <t  *> 

P  Pf 

when  i _ la  tha  ooaponant  of  1  in  tha  radial  direction.  Thua  wa  naad 

pr  p 

only  calculate  tha  dlaplacaaanta  for  a  ona-dlaana lonal  problea.  Wa  obtain 
valuea  for  v  and  <lpp*>  by  alaulating  tha  intaraetlona  between  a  partlcla 
and  a  large  nuabar  of  randoaly*or lanead  ona-dlaanalonal  addlaa. 


for  a  atralght  olroular  pipe  tha  addlaa  were  found  by  Hutohlnaon  to  ba 
approxlMtely  uni  fora  in  also  and  valoolty,  and  (Ivan  by 


1#  •  0.22  K 


and  •  Ut 


(3) 


where  la  tha  char  an  tar  latle  friction  velocity  and  f  is  a  single-phase 
friction  factor.  Tha  addy  llfatlaa  waa  alao  aasuaad  conatant  and  glvan  b> 


t#  -  1.6  i#/U# 


•tinea  tha  partlolaa  ara  aaall  thalr  aaan  valoolty  will  be  cloaa  to  tha  main 
gaa  velocity  and  ao  there  will  be  negligible  alip  in  the  axial  direction. 
Tha  interaction  tine  for  aaoh  addy  la  therefore  glvan  oy 


T  -  aln  <T#,  T^) 


where  la  tha  tlM  taken  for  tha  particle  to  croaa  the  eddy  In  the  radial 
direction. 

2.2  lauatlona  of  notion  for  aaall  oartlclaa 

To  oaloulata  tha  aaan  square  dlaplacaMnt  wa  oust  integrate  the  equation  or 
action  of  the  pertlale,  over  the  time  T  given  '.y  equation  (»>.  For  smell 
partlolaa  gravity  may  ba  neglected  and  tha  aqua  t  Ion  of  not  Ion  Is 


dU 


id 


d?  '  CD  H  •’g-^  'He  * 


)  M!n  '  U. 


for  aaall  partlclan,  the  Blip  in  the  axial  direction  is  usu.il:,  sir- 

aaalltr  t h k n  t  h  •  tddy  velocity.  For  •xanpl^,  f on  u  _  •  *  t  n  a r <\ 

G 

h  •  0.165  ■.  •  typical  eddy  valoolty  is  aeout  0.?->  m/s.  The  terminal 
velocity  of  a  partlcla  with  density  1 600  kg/is’  and  dlamet»r  .'■O  „r.  is  irour 
0.02  a/a.  Thua,  for  aaall  particles,  we  can  write 


■  dUp*  *v 

«-  *C0*»,T  <"o.  *  Vl°*  ‘°p«l  (6> 

for  the  velocity  oonponont  In  the  radial  direction.  Th«  naan  gaa  flow  hat 
M  horlgontal  ooaponent  and  no  UC|  •  (-1)?U#  where  q  la  a  rartdoa  nuaber. 
Hutch  I  neon  uaad  the  following  aquation  for  tna  drag  ooafficlant 
mc 

C0  -  jy  whore  C  •  0.116  (log,,  ROp)*  *  0.05»»  log,,  Rep  •  l.eej. 

P 

equation  (6)  waa  tntagratad  by  Hutch lnaon  at  al  to  give  tna  particle  velo¬ 
city  and  dlaplaoaaant  aftor  each  interaction,  and  hence  the  aean  aquara 
dlaplaoraent  par  lntaraotlon. 


2.3  tduatlona  of  notion  for  large  aartlclaa 

for  largo  partlolaa  tha  allp  In  tha  ailal  direction  will  usually  be  auch 
greater  than  in  tha  radial  dlreotlon.  for  eraapla,  tha  terminal  velocity 
of  I  I  ■  particle  with  a  danalty  of  1600  kg/«*  la  6  n/s.  Equation  (51 
thua  bacoaaa 


OIL  »««' 

■-i?*coS-s  ^-^p’ 4U. 


where  AU^  la  tha  velocity  dlffaranca  In  tha  axial  direction,  aaauswd 
constant. 

for  tha  radial  component,  th la  glvaa 
dU  ed  ' 

■if5  *°p,)‘ux 


Integrating  over  t laa  T  glvee 


«d  * 

-  >  log((-nq  U#  *  «p,)  -  Cp*  o-g  AUt  T 

Rearranging  and  using  d Inane  ion leae  vet  tables  V 
1 1  •  T/(wu#)f  s  -  •  1.  h  •  tpr/i#,  gives 


¥  -  S 


(S 


V 

o 
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3 


Miere  h  •  X 


Vp 


Intagrating  again  glvaa 
(S-Vj 


<S-Vj  -t,*4U  /U 

*  •  -  nu>  (’  -  •  1  1  #) 


(8) 


a  a 


Bacauaa  the  allp  la  larga  in  tha  axial  direction,  tha  lntaractlon  tlaa  win 
alwajra  be  tha  tlaa  takan  to  oroaa  tha  eddy  In  that  direction 

U, 

la  t  •  •  Thua  *i  “  JJJ  •  Thus  aquation  (?)  becomes 


S  -  (S  -  V  )  a 
o 


-A 


(9) 


and  aquation  (8)  baooaaa 


U  (S  -  (S  -  V  )  (1  -  a'*) 

»  •  S3,  -T1 


If  A  <<  ‘ ,  than  1 

h  - 

AUl 


a  *  •  A,  ao 


(10) 


(11) 


2. »  blffualon  ooafflolant  for  lar«a  garticlaa 
Equation  (2)  My  ba  wrlttan 
K,  .  .  2  <h*  >  1  ' 


■  .  " 

'51.  i  I  V 

*  i*i  1 


whara  tha  hj's  ara  tha  dlMnalonleas  displacements  calculated  from  tne  T 
alaulatlon.  It  has  alraady  baan  polntad  out  that  this  equation  applies 
only  whan  Individual  displacements  are  uncorralated.  For  large  particles 
this  la  unllkaly  to  ba  tha  case  aince  a  particle  moving  ir,  a  particular 
dlractlon  will  have  too  great  a  momentum  to  be  Immediately  llvertei  to  arc, 
other  direction.  In  general 


-  u 


K  -  |  1,'  i  (  I  \)' 
*  *  "  1-1  * 


,  »  * 

<«;•>  f  I  I  hi  hj 

n  "  1-1  j-i  J 

Ko  1  *  * 

<*•>  1  ii  jli hl  hl’J 


where  the  witlw  over  J  Is  understood  to  "wrap  around"  ths  value  of  l. 
K  .  M  M 

*  <iJ«>  if  E  E  <hi  hi-j> 

<n  *  "  l-i  j-i  1  J 


Thus  K 


K  II  II 

5°  I  I  e”*^”1  *  mint  equation  (A5) 

1-1  J-1 


Ko  I  * 
J-1 


Convert ln«  the  eu east  Ion  to  an  integral. 


K  -  Ko  /  e"*^’1  }dj 


I*  H 


-(M-1  )*, 


(12) 


If  h  Is  sufflolently  large,  this  gives  K  -  l^/A.  Otherwise,  aquation  (12) 

Z  .4U 

■ust  be  evaluated,  with  M  given  by  M  •  v — n-  .  The  value  or  Kq  may  be 

e  pi 

obtained  froa  the  10  slaulatlon  discussed  In  Section  2.1.  However, 
equation  (At)  shows  that  It  can  be  approx  lasted  by 


K 


o 


A  U. 

2 


.2 


4U*  *.  Ue*' 

Putting  v  •  ? —  gives  Kq  -  — 

e  z 


(  13' 
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thus,  for  a  sufficiently  long  pip* 


».  V 


3.  TWO  D1KM8I0MAL  SIHULATIOW  Of  PARTICLE  TRAJECTORIES 


To  Investigate  tho  validity  of  tha  above  aethesatlcal  troataent  of  the 
drop la t  aotlon,  a  taohnlqua  baa  baan  davtaad  In  which  the  particle  distri¬ 
bution  la  calculated  by  plotting  tha  trajectories  of  a  large  number  of 
partlolaa  aa  they  lntaraot  with  a  aucoaaalon  of  addlaa  of  specified  size 
and  veloolty. 


Tha  taohnlqua  la  described  In  aore  detail  by  Janes  at  al  (1980)  who  applied 
It  to  the  notion  of  snail  droplets  flowing  oo-ourrantly  with  a  gas  at>-eaB. 
Thay  uaad  tha  sans  description  of  tho  addlaa  as  In  section  2.1  or  this 
report  out  thay  treated  the  addlaa  as  two-dimensional  and  used  a  random 
nua bar  to  describe  the  orientation  of  eaoh  eddy.  They  used  the  following 
approxlnatlon  for  the  drag  ooefflolenti 


21/ 


He. 


0.4* 


because  using  this  equation  it  Is  possible  to  resolve  the  drag  law  into 
ooaponents  and  Integrate  over  the  Interaction  time. 


Boysan  at  al  (1982)  uaed  a  similar  technique  for  the  motion  or  droplets  in 
a  oyolone.  However,  Instead  of  assualng  constant  eddy  properties,  tney 
calculated  looal  values  fro s  the  following  equations 


0.3  k^^/c 


T#  •  0.3 1  k/c 


-  b  - 


Mwrt  k  and  c  ere  the  turbulent  kln«tlo  energy  end  dissipation  rats  pre- 
dlotad  by  a  k-c  model  for  the  single  phase  now,  and  q, ,  q,  art  pseudo* 
random  number  c  with  a  noraal  distribution.  The  drag  coefficient  waa 
asauaad  constant  during  any  interaction. 

The  techniques  of  Jaass  and  Boysan  aay  both  ba  applied  to  countercurrent 
flow  provided  the  equations  describing  the  partlole-eddy  interactions  are 
aodifled  as  shown  in  section  2.3  to  allow  for  the  large  vertical  slip. 
However,  for  partloleo  f’owing  lr.  s  straight  pipe  (not  close  to  the  wail) 
the  two  approaches  gave  very  stellar  results. 

S.  COMPARISON  Of  PBEDICTIOHS  Of  DlffUSIOH  MODEL  AMD  TRAJECTORY 
SIMULATION 

The  diffusion  Model  of  Seotlon  2  and  trajectory  alaulatlon  of  Section  3 
have  been  used  to  oaloulate  the  distribution  of  particles  dropped  verti¬ 
cally  on  the  centre  line  of  a  0.33  ■  disaster  vertical  column.  The  calcu¬ 
lations  were  oarried  out  for  partioles  of  disaster  2,«  an  and  1  mo,  density 
1600  kg/a'  and  gas  velooitlea  of  6  a/s  and  3  a/s 

Provided  the  particles  do  not  reach  the  walla  then  the  solution  to  the 
diffusion  equation,  equation  (1),  for  these  conditions.  Is  a  two- 
dlaenslonal  Oausslan. 

W(x.y)  -  jjj.  •  ’-<«• 

where  s'  •  2Kt 


Integrating  over  all  values  of  y  gives 


J _  #-x*/2o' 


which  is  a  Gaussian  with  rna  displacement  0. 


figure  1  shows  the  particle  distribution  at  z  -  3.6  a  predicted  ny  tne 
simulation  model  using  500  particles.  Clearly  this  is  also  Gaussian,  with 
0  given  by  the  ras  value  of  x^.  Thus  an  effective  diffusion  coefficient  k 
aay  be  obtained  froa  the  calculated  value  of  0  for  .his  distribution. 


f  lgure  2  shows  average  valuas  of  the  oorralatton  terms  <  hjh  >  calou.ated 
by  the  alaulatlon  model.  The  magnitude  of  these  taraa  decays  very  jiowiy, 
with  several  bund  rad  Interactions  required  bafora  the  oorralation  decays  to 
a  negligible  value.  The  predictions  wt  aquation  (16)  ara  also  shown  and 
agrea  qulta  wail  with  tha  alaulatlon  aodal. 

Consider  now  th'  diffusion  ooaffiolanta  oaloulatad  with  and  without  tha 

oorralation  tcraa.  In  Tabla  1  K  la  tha  diffusion  coafrtolant  obtainad 

o 

whan  oorralation  taraa  ara  nogleoted,  and  K  Is  tha  valua  obtainad  whan  thay 
ara  lnoludad.  for  tha  conditions  uaad  hara,  K  is  two  ordars  of  magnitude 
greater  than  KQ.  Tabla  1  alao  shows  that  tha  values  of  K  and  Kq  pradictad 
by  tha  simplified  equations  (17  and  18)  agree  with  valuas  obtained  from 
tha  10  and  20  simulations,  aapaolaily  for  the  larger  particles  (which  nsva 
a  higher  valua  of  A).  Hara  M,  tha  total  nuabar  of  interactions,  •  1 30  and 
130  for  tha  2. A  am  and  1.0  am  partlolaa  raspaotivaly. 

5.  C0HCLU3I0H8 

Tha  diffusion  aodal  of  Hutchinson  at  al  (1971)  for  snail  particles  in  a 
turbulant  gas  now  has  bean  extended  to  large  particles  in  a  vertical  pipe, 
where  tha  slip  between  tha  particles  and  tha  assn  gas  flow  is  large. 

It  was  shown  that  under  these  olrcuaatancea  the  displacements  resulting 
from  successive  pertlole-eddy  Interactions  are  correlated  ie  there  Is  a 
high  probability  that  suceesslvs  displacements  will  be  In  the  same  direc¬ 
tion.  These  correlation  effects  osn  result  In  diffusion  coefficients  which 
are  two  o-ders  of  magnitude  h.gner  than  would  be  expected. 


le  small  (la  A<<1)  than  tha  aquations  describing  tne  particle  eddy  .  '*r- 
aotlona  can  be  greatly  slmpllflod  and  tne  following  equation  was  dr-ived 
for  the  diffusion  coefficient.  Including  the  correlation  errect  just 

daacr  lbe-i. 
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l.  V 

TH7 


1  -  • 


-(■-1 )A  , 


M »»«.*•  M 


3. 


Particle  distribution*  calculated  with  this  diffusion  coefficient  agreed 
well  with  the  distribution*  predicted  by  *  trajectory  simulation,  Tor  a 
particular  set  of  oondltlons.  This  shows  that  the  mathematics  or  the  pro* 
cess  have  been  correctly  described.  However,  alnc>.  the  physics  or  Lhe 
partlole-eddy  Interactions  is  the  same  in  both  models,  this  aspect  has  not 
yet  been  tested. 


It  Is  intended  to  carry  out  experiments  In  a  perspex  test  section  using 
photographic  methods  for  determining  particle  position.  A  parametric  study 
will  also  be  carried  out  using  the  model  to  determine  the  limits  of  vali¬ 
dity  of  the  model. 
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WMMX  -  Calculation  of  oorralatlon  taraa 


equation  (9)  My  M  written  a  a  a  recurrence  ralatlon 


*1  -Sl’(Sl  -»!!)• 


"A 


(AO 


Applying  tMs  ralatlon  n  tlaaa  gives 

»n  -  (i  -  a"*)  I  s1.’t"’1>A  ♦  a "nAV<> 

If  n  la  aufflolantly  larga  tha  tara  e*^  ’/  aay  be  negloeted. 
•.^plying  aquation  (A1 )  “jo  i»n,  and  1-n-r, 


(A2) 


V  •  s„  <<-•”*)  *  V„  ,  e“A 
n  n  n-1 


“  Vr  '  Vr  l,V‘>  *  W"* 


Multiplying  and  auaatng  ovar  n, 

H  N  .  .  .  N 

l  Vn-r  -  I  •  VlVn-*-1*  •  O  *•  )  I  I  ,S  „ 

n-i  n  n-r  nil  n  1  n~r~'  ni,  "*>  n-r 

N  N 

♦  O-a’V  I  3nSn.r  ♦  e"A(i-e*A)  J  S  v 
n-1  n*i 

Tha  last  two  auaaatlona  aay  ba  naglectad  alnoa  the  S  ’s  are  uncorrelated 

n 

and  there  la  no  rorrelatlon  between  Sn  and  Vn_r_t ,  provided  r  >  o.  Thus 

N  N 

I  vvo-e'2A:  ;  ’*  (i-,*A)  l  vn  ,  s„  „ 
n-1  n  *  n-1  n  1  n  r 

-a  -»  N  n-1 

•  e  (1-e  *)  l  S  l  O-e  *)  S  a' 
n-i  n  r  l-i  1 


-A  ,  .  -(  n-i -l  )A 


-*-*  H  -(r-OA 


using  equation  (A2) 


-A, .  -A ,«  5 

a  (1-e  )  J  • 

n-i 


alnoe  the  S_'s  are  uncorrelated 
n 
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k  oonstant  defined  In  equation  (7) 

Cp  dreg  ooefflolent 

dp  pertlele  disaster 

n  dlaenslonlese  partlole  dleplaoeeert 

k  turbulent  klnetlo  erergy 

K  diffusion  ooe  flolent 

K,  diffusion  ooefflolent,  neglecting  correlation  terma 

1#  eddy  lengthscale 

lp  particle  displacement 

lpp  partlole  dlsplaoeaent  In  plane  perpendicular  to  gaa  flow 

a  partlole  eaaa 

N  total  nuaber  of  parttola-eddy  Interactions 

q  randoa  nuaber 

qia  q,  noraally  distributed  paeudo-randoa  nuaber  a 

R  pipe  radlua 

RCp  particle  Reynolds  nuaber 

S  dlaenslonlesa  eddy  veloolty 

T  Interaction  tlae 

T'  eddy  llfetiae 

U#  eddy  veloolty 

UQ  aean  gaa  veloolty 

Up  partlole  veloolty 

dlaenslonlesa  particle  veloolty 

V,  dlaenslonlesa  partlole  veloolty  at  start  or  interaction 

W  particle  density  distribution 

x,  y  co-ordinates  In  plane  perpend loular  to  tube  axis. 

Z  length  of  tube 

AUz  particle  slip  veloolty  in  axial  direction 

c  turbulent  dissipation  rate 

v  interaction  frequency 

gas  density 

e  atandard  deviation  of  partlole  distribution 

i(  dlaenslonlesa  interaction  tlae 

Subeorlpta 

x,  y  directions  perpend loular  to  tube  axis 

l,J,n,r  ith,  Jth,  nth,  rth  interaction 
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TABLE  1 


COHPMISOM  Of  DimiSIOM  COCrflCIOtTS 


%  <-> 

00  (./.) 

Ad^a/a) 

"•p 

CD 

* 

N  (for  Z  -  3.6  ■) 

K,  (aa*/s)  pradlctad 
by  ID  alaulatlon 

K,  pradlctad  by  aquation  03) 

K  (aa*/a)  pradlotad  by 
aquations  (12)  ant)  0  3) 

K  avaluatad  froa  pradlotlona  of 
2D  alaulatlon  audal  at  Z  *  3-6  a. 


2.* 

6 

9.5 

1500 

o.m 

0.005 

1 80 

0.37 

0.4? 

50 

43 
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1.0 

3 

5.9 

390 

0.59 

0.0174 

'30 

0.57 

0.73 

38 

4*) 


END 

DATE 

GtLMED 


